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ABSTRACT 

The passage of a weak vorticity disturbance through a reactive shock wave, or deto- 
nation, is examined by means of a linearized treatment. Of special interest is the effect, 
of chemical heat release on the amplification of vorticity in particular, and on the dis- 
turbance pattern generated downstream of the detonation in general. It is found that the 
effect of exothermicity is to amplify the refracted waves. The manner in which the imposed 
disturbance alters the structure of the detonation itself is also discussed. 
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INTRODUCTION 


When a pattern of vortieity in an otherwise uniform stream passes through a plane 
shock, it undergoes refraction and amplification, with the simultaneous generation of acous- 
tic and entropy signals behind the shock. For small amplitudes of incident disturbance, 
a linearized analysis is possible, and has indeed been carried out by several authors; see, 
for example, Ribner [1], McKenzie and Westphal [2] and Anyiwo and Bushnell [3]. The 
predictions of the linear theory and departures therefrom have been examined by Zang, 
Hussaini and Bushnell [4] via two-dimensional numerical simulations. 

The above studies have all been confined to chemically inert flows. With renewed 
interest in hypersonic vehicles, there has emerged a need for improved understanding of 
shock-turbulence interactions in reacting flows. This paper addresses a specific reactive- 
flow configuration, namely, a planar standing detonation wave (reacting shock wave). Such 
waves, provided they can be stabilized (still a matter under exploration), have been pro- 
posed as an alternative to the SCR AM JET concept in high-speed propulsion (O’Brien 
and K obayashi [5]). Our aim is to examine the interactions of the detonation with a 
small-amplitude shear wave as the latter is convected through the detonation. This is 
accomplished by means of a linear analysis, which can be thought of as an extension of 
Ribner’s [1] for an inert shock. Unlike Ribner, however, we are interested not only in the 
effect of the detonation on the disturbance, but also on the manner in which the shape 
and structure of the detonation are altered by the disturbance. 

THE GOVERNING EQUATIONS 

Figure 1 illustrates a stationary, planar, oblique detonation wave, into which flows an 
unburnt mixture of reactants from the left, with burnt products emerging on the right. The 
undisturbed flow is assumed to be uniform on either side, with the normal component of 
velocity supersonic ahead of the wave and subsonic behind. In the (.r, y) coordinate frame 
the detonation shock is represented by x = 0, the .r — axis being normal to the shock. (The 
analysis assumes that there are no variations in the z — direction.) The (£,?/) coordinate 
frame is aligned with the downstream flow, and the (£',?/') frame with the upstream flow. 

Let us now suppose that a planar, steady, shear disturbance, whose amplitude is small 
and depends only on the cross-stream coordinate 7/', is superimposed on the incident stream 
(Figure 2). The intent is to compute the perturbed flowfield behind the detonation and 
the perturbation of the detonation structure. 

It is convenient to nondimensionalize the system. Accordingly, the pressure, temper- 
ature and density of the gas, and the mass fraction of the reactant in the mixture, are 
referred to their respective values in the unrcactcd, unperturbed gas ahead of the detona- 
tion. Velocities are referred to the frozen sound speed in the unperturbed gas, and lengths 
to the characteristic wavelength of the disturbance. Then, the basic equations expressed 
in the downstream coordinates (£,?/) are 


(la) 


(P u )t + (P v )v = °> 


(ib) p( un t + vu v) + (ih)pt = °> 

(lc) p(uvt + vv n ) + (lh)Pt, = 0, 

'y — | 

( 1 J) p(uTt + wT,) — (»/>* + vp v ) = 


(le) l> = pT, 

and 

(1/) = ~ A > 


whore 


(^) 


AL 0 

0 

c 


K/7’ 


Horn, (u,v) is the gas velocity, V the reactant, mass fraction, and />, p and T the density, 
pressure and temperature respectively. A single irreversible reaction of type A -* D, 
governed by Arrhenius kinetics, has been postulated. (The details of the reaction scheme 
will only influence the effect of the disturbance on the detonation structure; the effect of the 
detonation on the disturbance depends only on the overall heat release.) The dimensionless 
parameters 7, E and n are, respectively, the ratio of specific heats, the activation energy 
and the heat release parameter. These are defined in terms of appropriate dimensional or 
reference quantities by the expressions 


7 = Cp/C vt o = QY/C P T, E = E/RT, 


where Y is the reactant mass fraction and T the temperature of the fresh mixture, while 
Q is the heat released when a unit, mass of the reactant is consumed. The dimensional 
activation energy is denoted by E , the gas constant by R and the specific heats at constant 
pressure and volume by C p and C,, respectively. Finally, Lq is the characteristic wave length 
of the disturbance, c the sound speed in the unlmrnt gas, and A the pre-exponential factor 
in the Arrhenius rate law. 

Typically L 0 is much larger than the thickness of the detonation, so that on the scale 
of the disturbance, the detonation can be treated as a discontinuity in an inert flow. When 
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examining the effect of the detonation on the disturbance, therefore, equation (If) can be 
replaced by 


Y = 1 ahead of the detonation, 0 behind the detonation. 

However, (If) will need to be reinstated, and the entire equation set ( 1 ) rescaled with 
respect, to the detonation thickness, when assessing the effect of the disturbance on the 
detonation structure. 


EFFECT OF DETONATION ON DISTURBANCE 


With the detonation treated as a discontinuity, the generalized Rankine-Hugoniot 
conditions yield the following expressions for the state of the gas behind it: 


(2a) 


1 + 7 A / 2 
(7 + 1 )M 


1 

(7 + 1) 


'(M 2 - l ) 2 

M 2 


1 1/2 


2(7 + l)o 


(2b) 


p = M/m , 


(2c) 


T — 1 -f rv + 


MlW'-™')! 


(2rf) 


P = pT 


M 


m 


1 + fir + 




(2c) 


to = IF or equivalently, m tan <p = M tan#. 


Here, m is the normal velocity and w the tangential velocity behind the discontinuity, while 
M and IF are the corresponding quantities ahead. Observe that M is also the normal 
Mach number in the unburnt region. The minimum value Mcj of M, corresponding to a 
Chapman- Jouguet wave, is given by 

(3) M C j = [l + (7 + l)n + [{1 + (7 t 1>>} 2 - 1] ,/2 ] 1/2 . 

A graph of Mcj against o is shown in Figure 3. 

Let the perturbed upstream flow be characterized by 

(4 a) U = U 0 + Ut, Ui«Ua, 

and correspondingly (see Figure 4), 

(4ft) # = M = M q + M] , W = W 0 + IF, , 
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wliiln all other state variables upstream remain unchanged. (The precise form of U j will 
he specified later.) It is now a simple matter to linearize the jump- relations (2) to obtain 
the following expressions for the disturbances immediately behind the detonation front 


(identified by the suffix /): 

(5a) 

u xf /via = —dcos<po[(Ut /Un) — (0| /?•) tan </7 0 ] 
+ sin </?o [^1 »' + (Ui/Uo) fallen], 

( 5 b) 

Vif/m 0 = d sin </>o[(Ui fUo) — (#1 /r) tan y> (J ] — 6 X sec 
+ cosy>o[^i r + (U\/Uq) tanv9o]> 

(5c) 

pi f = 7 [2 - (1 - d)/r][(Ui/U 0 ) - (^,/r) tan^o], 

(5c/) 

T xf = (7- 1)M, 2 (1 + d/r 2 )[(U| /U o )-( 0 i /r) tan y> 0 ]. 


The constant r is simply the ratio of the normal speeds across the detonation, 
(G) r = Mo/»»o, 


while d is defined by 


(7) 


d = 


(7 + 1 )M 0 vt 0 


1 - 7 K + 


K 


{(Mo - l) 2 -2(7+ 1M/ 2 } 1/2 J 


The expressions (5) act as boundary conditions for the perturbed flow downstream 
of the detonation. The governing equations for this flow are simply the Euler equations 
(la-e). Upon linearization about the basic downstream state these reduce to 

(8) (l - + Vv/ = 

where is the perturbation stream function, /< the Mach number of the undisturbed 
downstream flow, defined by 


(9) /< = (>'S/r„) ,/2 , 

and fii the vorticity behind the wave. It can be shown that I2| depends on the cross-stream 
coordinate 7/ alone, and can therefore be evaluated at the location of the undisturbed front, 
again identified by the suffix /, i.e., 


fl. 


1 

P\ 7 1 


+ u 1 r/ 


J/ 


( 10 ) 
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Once 0 is known, the velocity components are given by the expressions 


( 11 ) 

«i=Vv ”1 = -(1 - /*>/’( 

while pressure, 

temperature and density can be found from 

(12 a) 

Pi = Ptf + 7/»o«o( M i/ - «i), 

(12 b) 

r, = t„ + _ Pl/ ), 


7/’o 

(12 c) 

P\ = (l/7o)(pi — poTj). 


Now, depending upon the downstream Mach number //., equation ( 8 ) is either elliptic 
(// < 1 ) or hyperbolic (/*>!). In turn, ft is an increasing function of the inclination 6 0 of 
the incident stream for a given upstream normal Mach number M 0 . In fact, there exists a 
critical value 9 C of 0 O , given by 

(13) tan0 c = [{2(1 + o) + (7 — l)Mfi — (7 + l)ntp}/( 2 M^)] ^ , 

such that fi > 1 for 9 > 9 r and /i < 1 for 9 < 9 r . Figure 5 shows a plot of 9 C as a function 
of the upstream normal Mach number A/o (more specifically, Mq /Mcj') for various values 
of the heat-release parameter a . Observe that for a CJ — detonation, the downstream flow 
is supersonic for all 9 0 > 0 . 

Thus far, the formulation has been quite independent of the functional form of the 
incident disturbance. Following Ribner [ 1 ], let us take a sinusoidal perturbation, i.e., 

(14) U\ - fUo cos A'V, 

where e characterizes the perturbation amplitude and A-' is the incident wave number. The 
downstream wave number k, given by 

(15) k’ cos 9 0 = k cos if 0 , 
must also satisfy the matching requirement 

(16) k 1 if = ki) 

at the front. Therefore, 

(17) U 1 = fUo cos An/ at the front. 

The corresponding perturbation in the front inclination may be taken to be 

(18) 9 1 = f(«cos An; + Asin An/), 


C 


where a and b are to be determined. 

The problem now reduces to the following: using the prescriptions (17) and (18) in 
the boundary conditions (5a, 6), solve (8) for V’’ Then, uj and V\ can be found from (11) 
and the remaining variables from (12) upon using the conditions (5). Only the results 
are given below since computational details are similar to those in Ribner [1]. (We have 
elected to display the velocity and pressure fields only; temperature and density follow 
from {12b, c) in a straightforward manner.) 


For // < 1 the solution is 


(19a) 


ui /(eU 0 ) = S cos[J-j,(?y — x tany>o) + 

+ P{r)cos[k y (y - xtan^) + 6 r ], 


( 106 ) 


Vi/{eU 0 ) - ftP(x)sm[k y (y - x tan^J,) + 6,], 


(19c) 

where 

( 20 ) 


pi /( eU 0 ) = -^MoP{x)sccApn cos[fc 9 (y - xtanepo) + b r ], 


ft = (|1 - l< 2 1)' /2 , ftl = 1 -/' 2 cos 2 <po, 


( 21 ) 

and 

(22a) 


k v = k cor 0 ( >, 


P(x) = [1 /(ftr)J(Cf + C| ) 1 exp( — k y ftx / ft „) cos 0 O , 


(226) 


S = {A 2 + Z? 2 ) 1 ^ 2 ( 1 /r ) cos 8 $ , 


(22c) 


/>, = tan l ( — B/A), 


(22 d) 


— tan 1 [(C'i ft — C 2 tan ^p^)/{Cift + Ci tan <po )]> 


(22c) 


<fo = - tan 1 [( 11 /ft) 2 cos 2 <p 0 tan <p 0 ]- 


The constants A,D,C\ and C 2 appearing above have been defined in the Appendix, along 
with the constants a and b introduced in (18). 
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For //. > 1, 

«,/(€l/o) = -S’ ros[fr v (j/ - ar tan V?o)] 

+ F cos[ty(j/ - .T tan (^o)], 

v\/(fU 0 ) = ftPcns[k y (y — x tan <^o)], 

Pi/((Uo) = sec <p a cos[ky(y - x tan^)], 


P — (C*j/r)(cos t?o sincj)/cosy3o, 
w = cot, -1 /?„, 

5 = (A/r)cos0 o , 

An examination of (19) and (23) reveals that the velocity field is made up of a shear 
component of amplitude S and a potential-flow component of amplitude P. The shear 
component is a plane wave which carries all the vorticity and propagates in the stroamwi.se 
direction. The potential component, as Ribner [1] has pointed out, is acoustic in character 
and is the sole contributor to the pressure disturbance. It also travels as a wave, but in a 
direction inclined at an angle tp ' to the normal, with amplitude remaining constant for the 
supersonic case and attenuating exponentially for the subsonic. The subsonic flow involves 
phase shifts 6 a and 8 p , while the supersonic flow does not. 

The amplitude 5 and the phase shift, are plotted in Figure G(a,b), against the wave 
inclination 8q, at Mq/Mcj = 1.5 and for various values of the heat release o. The most 
dramatic effect occurs at the sonic point,, 6 — 0 C , where the amplitude increases sharply. 
Figure G(c) displays the maximum amplitude as a function of at for four different values of 
the incident Mach number, and reveals that for rv beyond 0.5, the graphs are essentially 
linear. In Figure 7(a,b,c) the corresponding graphs for P , the amplitude of the potential 
component (evaluated at the shock for subsonic flow), and the phase shift 8 p are shown. 
The behavior is similar to that in Figure G, with one exception: whereas the maximum 
amplitude of the shear component increases with increase in Mq/Mcj (Figure Gc), that of 
the potential component decreases (Figure 7c). 

Because of the nonuniform pressure jump across the front, it develops ripples, and 
the resulting deflection of the front from the undisturbed position x = 0 (obtained by 
integrating the front-deflection angle 0] ) is found to be 

(25) Ay = [,l/„/(VA/ 0 )l(„ 2 + I, 1 ) 1 ' 2 (■<*(»■,„ + 


(23„) 

(23b) 

(23c) 

where 

(24a) 

(246) 

(24c) 

(24c/) 
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whore k' is the wavenumber of the incident, disturbance and «, b (introduced in (18)) are 
defined in the Appendix. The quantities ( a 2 + 6 2 ) ,/2 and S tlonl have been plotted in 
Figure 8(a,b) against the wave inclination 0 o for different values of M 0 /M C J • Again, 
heat release makes itself felt most significantly at the sonic point. In Figure 8(c), t,he 
maximum of (a 2 + 6 2 )'/ 2 is plotted as a function of o for various Mach numbers. Unlike 
the corresponding graphs of Figures 6(c) and 7(c), the variation with Mach number is 
minimal here. 


EFFECT OF DISTURBANCE ON DETONATION 

Let us now examine the alteration in the structure of the detonation front caused 
by the passage of the disturbance. This will be done in the asymptotic limit of large 
activation energy. First, it is useful to recall the structure of the undisturbed oblique 
detonation, displayed schematically in Fig. 9. The structure consists of a lead shock S, 
across which the state of the mixture suffers a jump but the mass fraction retains the 
upstream value of unity since the shock is inert. Behind the shock lies the induction zone 
I in which the reaction is weak, and therefore, the state is undisturbed to leading order. 
The induction zone terminates in the fire zone F, which is a thin region of intense reaction. 
Behind the fire zone the mixture is fully burnt. The broadest feature of the detonation 
structure is the induction zone, whose thickness L/ effectively determines the detonation 
thickness. On the Lj scale, the shock is a discontinuity of course, but so is the fire zone in 
the limit E -+ oo. It can be shown that L, has the asymptotic form 


(2G) 


Li 


_ 5^ e K/7’o 

AE 


where To is the undisturbed temperature immediately behind the lead shock. In the large 
E limit, it is possible to give an analytical treatment of the induction zone, in which 
state variables change by 0(E ~ } ) and the position of the fire zone is characterized by a 
singularity in the induction-zone solution; the so-called thermal runaway. We now show 
that the primary effect of the applied disturbance is to produce relatively large ripples in 
the position of the fire zone. 

In the induction zone, it is convenient to adopt the coordinate system (£,fj), aligned 
with the undisturbed flow immediately behind the lead shock (Fig. 9). The dependent 
variables are denoted by hatted quantities as well, and the subscript zero corresponds to 
undisturbed flow. Also, the induction-zone thickness Lj is chosen as the reference length. 
Then, equations (la-f) apply again, provided all variables are hatted. Also, the reaction 
rate A now has the form 

A = ^±Y,-rn\ 
c 

which, in view of (2G), can hr cxprrsrd as 


(27) 



E 

T. ' 
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As before, the state of the gas behind the lead shock can be computed by means of the 
jump conditions (2), provided a is set to zero there. 

Although we have chosen to refer all lengths to Lj in the induction zone, it must be 
recognized that along the shock, i.e., in the y-direction, the proper length scale is still L 0 . 
Therefore, on the Lj scale, y-derivatives must be of order L//L 0 , and hence negligible to 
an excellent approximation. Thus, for any dependent variable F, 

F y - F^sin^o + Fycosifo = O(I//L 0 ) < 1, 

whence 

Fy ~ - tan tfoF^. 

The governing equations are thus rendered quasi-ono-dimensional in the induction zone. 

Let us now consider a distinguished limit in which the applied disturbance amplitude 
is T^jE. Then, for E — ► oo, the governing equations reduce to 


(28«) 

«o fy + - tan^oiip = (), 

(286) 

+ (1/7)7^ = 0, 

(28c) 

Po'iov'i - ( 1 / 7 ) tan <foP^ = 0, 

(28<7) 

- . rf* 7 - 1 . r 

PqUoTi - u q p~ = ac. , 

s y s 

(28e) 

p = f>[)T r + Top . 

In the above equations, primes denote perturbations. It is a simple matter to extract from 
these equations, by elimination, a single equation for the temperature perturbation T' , i.e., 

(29) 

= rP", 

where 


(30) 

T («7](l/T)-/< 2 n 
V ? 'oPo / 1-/1?, 

and is the normal Mach number of the undisturbed flow immediately behind the shock. 

We remark parenthetically that for T to be positive, /*„ < 1 /^ 7 , i.e., the normal flow speed 
immediately behind the shock must be less than the isothermal sound speed, a well-known 
restriction (see, for example, Clarke [G]). It is now a simple matter to integrate equation 
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(29) provided the value of the temperature perturbation at the shock, is known. 

Then, one obtains 

( 31 ) f' = -ln[e- f *-r£]. 


The actual profile of T' is of less interest than the position at which it becomes singular. 
This position is given by 

( 32 ) if = 

and defines the location of the fire zone. For the undisturbed detonation, T s == 0 and 
the fire-zone location is given by £ /0 = 1/r. The shift in location is therefore £/ - 1/r. 
Note that the shift is an 0(1) quant itity, on the scale of the induction length, while the 
disturbance amplitude causing the shift is 0(E '). The applied disturbance, therefore, 
has a pronounced effect on the internal structure of the detonation wave. This is illustrated 
in Figure 10 for a sinusoidal temperature perturbation at the shock. Once V is known, 
equations (28a-e) can be integrated, subject to appropriate conditions at the shock, thus 
yielding the profiles for all the variables in the induction zone. 

CONCLUSIONS 

On the whole, then, the effect of exothermicity is to amplify the refracted waves. This 
result, although derived here only for steady waves, contradicts that of Kumar, Bushnell 
and Hussaini [7]. However, they simulated numerically the interaction of a disturbance 
with an oblique shock with equilibrium chemistry , and found that chemistry led to a slight 
decrease in amplification. We believe that their result is biased by the manner in which 
they prespecify the degree of reactedness at each station. It would be more appropriate to 
compare our results with a direct numerical simulation of the reacting shock, undertaken, 
for example, by extending the approach of Zang, Hussaini and Bushnell [4] to reactive 
flows. 
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APPENDIX 


In this Appendix are listed the various constants introduced in the text. 


A = scctpo + 2 (r - 1) cos <p n + («/r)(r - l) 2 sin^o, 


D = {h/r)(r - l) 2 sin<^o, 


a 


{ 


..CE+UF 
1 17nr[yr 
r C' + GF' 
E'+air 


for // < 1 , 
for // > 1, 


/> = 


CF-DE 

7 C 2 + Z > 2 ’ 


C — \d + 2r 1 (r tan tpn — [(r — l) 2 //? 2 ] sin <£> 0 cos <£> 0 , 


C — 1 — d — 2[l -f (r — 1) cos 2 <£>o], 




D' = (r-l)[l + (r-l)oos 2 y, () ], 


E = 1 + rf + 2(r - l)(/? 2 //?, 2 t ) <-os 2 <^. 0> 


E' = (r - 1 ) 2 sin cos y?„ - (d -f r) tan v? 0 , 




F' = 2(r — 1) sin cos^> 0 , 
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_ 1-/7 tan y? 0 
/? + tan v?o ’ 

Ci = ar~ l D' — F ' , 

C 2 = tr-'D'. 
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Figure 1. Illustration of a stationary, planar, oblique detonation wave and the 
coordinate frame. 


L5 



Figure 2. Schematic of a planar, steady, shear disturbance superimposed on 
the incident stream. 
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Figure 3. Plot of the Chapman-Jouguet Mach number Mqj versus the heat 
release parameter a. 
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Figure 4. Schematic of the perturbed upstream flow showing the nomencla- 
ture. 
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Figure 5. Plot of the critical angle 0 C versus the upstream normal Mach 
number M 0 / M CJ for a = 0, 0.1, 0.5, 1, 5. 
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Figure 6b. Plot of the phase shift 5, versus the wave inclination angle 0 O , at 

M n / Mr] = 1.5. for a = °> 01 > °- 5 ’ 5 - 





Figure 7a. Plot of the amplitude P versus the wave inclination angle 0 O , at 
M 0 / M CJ = 1.5, for a = 0, 0.1, 0.5, 1, 5. 
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Figure 7c. Plot of the maximum amplitude of P versus the heat release 
parameter ot, for M o / Mqj = 1-5, 5, 10, 14, at 0,.. 


Figure 8a. Plot of V (n 2 + b 2 ) versus the wave inclination angle 0 O , at 
M 0 / M CJ = 1.5, for a = 0, 0.1, 0.5, 1, 5. 






Figure 8c. Plot of the maximum of V (a 2 + b 2 ) versus the heat release 
parameter a, for M 0 / M CJ = 1.5, 5, 10, 14, at 0 C . (Graphs for M 0 / M CJ of 
5, 10, and 15 are indistinguishable). 
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of the detonation itself is also discussed. 
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